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Abstract. This survey article deals with a delta potential - also known as a 
point scatterer - on flat 2D and 3D tori. We introduce the main conjectures 
regarding the spectral and wave function statistics of this model in the so- 
called weak and strong coupling regimes. We report on recent progress as well 
as a number of open problems in this field. 

1. Introduction 

The point scatterer on a torus is a popular model to study the transition between 
integrable and chaotic dynamics in quantum systems. It rose to prominence in 
the quantum chaos literature in a famous paper of Petr Seba [J5]. The model 
hrst appeared in solid state physics [TTJ in the 1930s to explain electronic band 
structure and conductivity in solid crystals. Many applications arose in nuclear 
physics throughout the 1960s and 1970s, see for instance [6]. The purpose of this 
article is to give an introduction to this important model which belongs to the class 
of pseudo-integrable systems and to report on some recent progress in this field. 
The reader will also be introduced to some important open problems. 

1.1. Kronig- Penney model. In 1931 Kronig and Penney [TT] studied the quan- 
tum mechanics of an electron in a periodic crystal lattice with the goal of under- 
standing the conductivity properties of solid crystals. 
They introduced the periodic ID Hamiltonian 

d 2 

(1.1) H KP = -— + VoY,X h zz ] (x-k), V Q >0, 0<a<l 

fcez 

where x denotes the characteristic function. According to Bloch theory, we have 
the decomposition 



L 2 (R) = f d9Ue, 9 £ [0,2tt) 



where "He is the space of quasiperiodic functions with quasimomentum e l6 \ 

He = {ip £ C°°(R) | ip(x) = e ifc V(x + k), k£ Z} 

Let us consider the special case of periodic boundary conditions 9 = 0. To simplify 
the Hamiltonian Hkp\-h it is convenient to take the limit 

a \ 0, a = Vqci = const. 

Let / £ C°°(R/Z). The calculation 

/a/2 
f(x)dx -> a/(0), a\0 
-a/2 
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shows that the Hamiltonian Hkp\u converges in the distributional sense to a 
singular rank-one perturbation of the ID Laplacian 

d 2 d 2 

(1.2) -—2 +VoX[_« ) |] H a =-^2 + <*(£<), •}$), a \ 0. 

The operator i? Q can be realised rigorously by using von Neumann's self-adjoint 
extension theory. We will be interested in studying the analogues of the operator 
H a on 2D and 3D tori. 

1.2. Seba billiard. Let R be a rectangle with side lengths a, b. We define the 
aspect ratio of R as the quotient (a/b) 2 . In a 1990 paper [19] Petr Seba studied the 
operator 

(1.3) H a = -A + a(S X0 ,-)S X0 , x £R,aeR\{0} 

on a rectangle with irrational aspect ratio and Dirichlet boundary conditions. Seba's 
motivation was to find a quantum system which displayed the features of quantised 
chaotic systems such as quantum ergodicity and level repulsion, yet whose classical 
dynamics was close to integrable. 

In the present paper we will deal with the closely related system of a point 
scatterer on a torus (which means periodic boundary conditions), however, the 
results which will be presented can probably be easily extended to rectangular 
domains with Dirichlet or Neumann boundary conditions. 

Acknowledgements: I would like to thank Zeev Rudnick and Stephane Non- 
nenmacher for many helpful comments and suggestions that have led to the im- 
provement of this paper. 

2. Quantisation of a point scatterer 

2.1. Self-adjoint extension theory. We consider a rectangle with side lengths 
2na, 2ir/a, where a > 0, and identify opposite sides to obtain the torus T 2 = 
M 2 /2tt£q where Cq = Z(a, 0) © Z(0, 1/a). We want to study the formal operator 

(2.1) H a = -A + aS Xo , ael\{0}, i eT 2 . 

To treat H a rigorously we will employ von Neumann's theory of self-adjoint 
extensions. For an introduction to this standard machinery see |27j . The main idea 
is to restrict H a to a domain where we understand how it acts - functions which 
vanish at the position of the scatterer and therefore do not "feel" its presence. 

We denote by 

(2.2) D Q = C%°(T 2 \{x Q }) 

the domain of C°°-functions which vanish in a neighbourhood of xq. Clearly 
Ho\d = — A|u . We denote — A = — A|.d . The restricted Laplacian — A 
is a symmetric operator, however it is not self-adjoint. By restricting —A to 
the domain Do we are enlarging the domain of its adjoint. Therefore we have 
Dq = Dom(— Ao) C Dom(— Aq). A simple computation of the adjoint — AJ shows 
that its domain is given by 

(2.3) Dom(~A*) = {/ G L 2 (T 2 ) \3CeC: Af + CS xn G L 2 (T 2 )}. 
We have the following definition. 
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Definition 2.1. The deficiency spaces of a symmetric densely defined operator A 
are given by the kernels 

(2.4) K± =kcr{^*±i}. 

The deficiency indices of A are defined as n + = dim/C + and n_ = dim/C_. If 
n + = ra_ = 0, then we say that A is essentially self-adjoint. 

For A ^ er(— A) denote by G\{x,y) the corresponding Green's function, namely 
the integral kernel of the resolvent 

(2-5) -^—f(x)=[G x (x,y)f(y)dy 
and therefore we have the following distributional identity 

(2.6) Gx{x , y ) = -±— Sy (x). 

Indeed, if we compute the deficiency elements of — Aq we have to solve 

(2.7) = (-A*±i)f = (-A±i)f + C5 Xo 

for some C e C. This shows that the deficiency spaces are spanned by the Green's 
functions G±i(x,Xo). We thus have 

(2.8) Dom(-A* )=Th® ± C{G i }®±£{G- i } 

where the orthogonal decomposition is with respect to the graph inner product 
(/, h) g = (/, h) + (Aq/, Aq/i) and the closure is taken with respect to the associated 
graph norm ||/|| g = (/,/) ? . 

The following theorem is due to von Neumann. 

Theorem 2.2. Let A be a densely defined symmetric operator. If A has deficiency 
indices n+ = n_ = n > 1, then there exists a family of self-adjoint extensions 
which is parametrised by U(n), the group of unitary maps on C". The domain of 
the extension Ajj is given by 

Dom(Au) = {fe Dom(A*) \f = g+ (G h v) + (G_ u Uv) , 

(2.9) 

g e Dom{A),v e C"} 

where G u G_i are the vectors whose entries are the deficiency elements and the 
closure of Dom(A) is taken with respect to the graph norm of A. The operator 
Au = A*\ Dom ( Au ) is essentially self-adjoint. 

In our case we have for ip e (— ir, tt] 

D v = Dom(-A v ) ={/ e Dom{-A* ) \ 

(2 ' 10) / = g + cGi + ce^G-i, g eTh,ceC} 

and — A v is the restriction of — Aq to this domain. Functions / e Dom{— A v ) 
satisfy 

3GeC: Af + CS X0 e L 2 (T 2 ) 
and near xo they have the asymptotic 

/(,) = G(^^ + tan|) +0( l). 

The extension associated with the choice tp = tt is just the self-adjoint Laplacian 
on G°°(T 2 ). We will be interested in studying the extensions for tp e (— tt,tt). 
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2.2. The quantisation condition. An orthonormal basis of eigenfunctions of the 
Laplacian on T 2 = M 2 /2ttCo is given by the complex exponentials 

Mx) = —e i <*>® where £e£ 
and C is the dual lattice of Co : 

(2.11) C = {£_ e M | V?/ € C : (£, 17) G Z} = {(m/a, rca) | m, n e Z} 
The eigenvalue of ip£ is given by |£| 2 . We introduce the set of distinct norms 

j\f = {0 = no < ni <n 2 < ■■■} 
and denote the multiplicity of n 6 Af by 

(2.12) r £ (n) = I n=|£| 2 }. 

If £ is a rational lattice, i. e. if a 4 <= Q, then the multiplicities can be large and we 
have the bound (see e.g. [TB], Lemma 7.2) 

(2.13) r £ (n) <C e n £ . 

We have the following lemma which can be found in the standard literature on 
point scatterers, or in the appendix to the paper [17] . 

Lemma 2.3. Let tp G (— tt,it). We have that X (fc o{— A) is an eigenvalue of —A v 
on T 2 iff 

(2.14) f>£K) (^ZTx ~ ^Tl) = c ° tan (^/2) 



wh 



ere 



Co 



E 
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and the corresponding eigenfunction is a multiple of the Green's function G\{x,xq) 
for which we have the L 2 -identity 

1 00 1 

(2.i5) GxiXtXo)= ^ E ei< ^ 0> - 

The eigenfunctions in the lemma are not the only eigenfunctions of —A v . If 
fc(nj) > 1 (which in our case happens for all rij > 0), then nj appears in the 
spectrum of — A v with multiplicity rc(rij) — 1. The associated eigenfunctions are 
superpositions of Laplacian eigenfunctions and vanish at xq an d therefore do not feel 
the effect of the scatterer. In the present article we will only be interested in the new 
eigenvalues, which are solutions to (|2.14[) . and the associated new eigenfunctions. 

3. Weak coupling vs. strong coupling 

Two different quantisations appear in the literature on point scatterers, and we 
will follow the terminology used by Shigehara et al. in the papers (2TJ [22j [23j [24] 
in referring to the two models as weak and strong coupling. The purpose of this 
section is to explain how these different quantisations arose in the literature and 
give a derivation of the strong coupling quantisation on a 2D torus. 

In his famous paper on wave chaos in a singular billiard [19) Seba considered a 
point scatterer on a rectangle with irrational aspect ratio and Dirichlet boundary 
conditions. Seba computed the spectrum of — A v by solving the equation (|2.14p 
numerically and a plot of the level spacing distribution seemed to suggest level 
repulsion. 

In 1994 Shigehara came to investigate this question and following careful numer- 
ical investigations he observed [21] that the eigenvalue spacing distribution of the 
self-adjoint extension — A v seemed to coincide with that of the Laplacian, which 
according to the conjecture of Berry and Tabor [T] is believed to be Poissonian. 
Shigehara observed that the apparent "weakness" of the point scatterer in Seba's 
quantisation could be corrected by adjusting the parameter ip in a suitable way as 
the eigenvalue A tends to infinity. 

The operator can be realised by employing the self-adjoint extension theory 
discussed in the previous section. The formal operators 

(3.1) -A + aS X0 , «eM\{0} 
are associated with the family of extensions 

(3.2) -A v , pe(-7r,7r). 

It is well known that in ID there is an exact relation which links the physical 
coupling constant a / and the parameter ip € (— 7r, 7r). The situation is more 
complicated in 2 or 3D. A relation can be derived (at a physical level of rigour) from 
the scattering problem for a spherical scatterer (cf. for instance [3D]) by shrinking 
its diameter to zero. 

One obtains a relation which in contrast to the ID case contains a logarithmic 
divergence in the spectral parameter A 

(3.3) - = d tan % + C 2 log A 

a 2 

for certain real constants C%, C2 which depend on the domain. 
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So a fixed choice of tp corresponds to weak coupling, 

1 



(3-4) a ~ 

C 2 log A 

On the other hand a fixed physical coupling constant a — const. ^ would require 
a renormalisation of the parameter ip, which means tp — tp\ should be allowed to 
depend on A as A — > oo, 

(3.5) Citan^ - -C 2 logA. 



The renormalisation condition (13.5[) is equivalent to a different quantisation con- 
dition for a point scatterer (cf. for instance [4]) which only takes into account the 
physically relevant energies in the summation, 



(3.6) 



|n 3 — \\<\ s \ J 



1 



— A n j + 1 J a 



where a =/= is the physical coupling constant. We require the following lemma 
which gives an asymptotic for the sum over the energies outside the interval [A — 
A 5 , A + A 5 ]. 

Lemma 3.1. Consider a general torus T 2 . Let S € (#,1), where 8 — We have 
the asymptotic 



(3.7) E ^^-^=-^ + 0(1) 

Proof. We have 



(3.8) '" ' 



E ^){-^--- 

n ? - — A n 7 



A ?i j + 1 



0(1). 



We will use the circle law 

(3.9) J2 r ^=™ + °^ e+£ ) 

rij <x 

where the best known exponent 8 = ^ is due to Huxley [9] and the optimal 
exponent is expected to be 8 = \ (Gauss's circle problem). 

Summation by parts allows us to compare a lattice sum with an integral 

(3.10) E /(Kl 2 ) = 7r f P f(t)dt + 0(a e f(a)+(3 s f(l3)) + o( f \.f(t)\t e dt\ 

where / is some differentiable function. 
We obtain for the first sum 

i n r r 1 n . 



(3.11) X+X s <rtj 1"J- A "jJ Jx+X 5 {X ~ A X \ 



(1 -<5)7rlogA + 0(A e -' 5 ). 
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Similarly for the second sum 

E i/ ^){ 5 -h-^}-r A, {^-i}' fa+o(i) 



(3.12 



nj <\-\ 6 



((5-2)7rlogA + 0(l). 



□ 



With the help of the Lemma above we can now show that the quantisation (|3.6j) 
is in fact equivalent to the strong coupling quantisation given by the renormalisation 
condition (13.51). 



We obtain 

°° f 1 1 1 

(3.13) g r ^)|__^j M]ogA + _ + 0(1) 

which forces the renormalisation 

(3.14) c tan ~ — rr log A. 

4. Statistics of wave functions 

Let {Aj} denote the solutions to the spectral equation (12. 14)) (where tp is fixed). 
The eigenvalues {Aj} interlace with the norms {rij} in the following way 

(4.1) Aq 5 < no < Af < m < ■ ■ ■ < Aj < nj < ■ ■ ■ 

We denote the solutions to the strong coupling quantisation condition Q3.6P by 
{A** r }. Similarly as above the eigenvalues {Aj* r } interlace with the norms {nj}. 
The corresponding eigcnfunctions are the Green's functions {G\st r . }. We are inter- 
ested in the statistical behaviour of the eigenfunctions {G x *} and {G\at r .} in the 
limit as the eigenvalue tends to infinity. 

Let a £ C°°(S*T 2 ) where S*T 2 ~ 't 2 x S 1 . Let (x,<j>) £ T 2 x S 1 and expand 
a £ C°°(T 2 x S 1 ) into a Fourier series 

(4.2) a(x,0)= £ 6(C,*)e i<C '" >+iW 

In order to quantise the classical symbol a it is convenient to realise the torus as 
the quotient T 2 = C/27r>Co[i], where 

C [i] = {xi + ix 2 | (xi,x 2 ) £ C }. 

For £ £ Co we define by £ the corresponding element of Co [i] . We associate with the 
symbol a the zeroth order pseudo-differential operator Op(a) : L 2 (T 2 ) ~ > L 2 (T 2 ) 
defined on the Fourier transform side by 

(4.3) (0K^/)(£)= ^ /(e-C), ^0 
and 

(4-4) (Opto/) (0) = £ ^' °)/(^ " 
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where 

(4.5) /(*) = £/(£)e i<£ ' x> - 

Let A € M \ <j{— A). Denote the L 2 -normalised Green's function by g x = 
Ga / II Ga II 2 - We would like to study the behaviour of the matrix elements 

Op{a)g x j , g x v 

in the limit as j — > oo. Similarly we would like to study the matrix elements of 
Op(a) in the strong coupling quantisation, 



(Op(a)g x .t r . , g x ,tr.) 



We first study the behaviour of the eigenfunction in position space. If we take a 
classical symbol on position space a E C°°(T 2 ), then the operator Op(a) is simply 
given by multiplication 



(4.6) (Op(a)g x ,g x ) = / a\g\\ dfi. 

Jt 2 

We have the following theorem, which is proven in |17j . In the paper the theorem 
is stated for the specific case of the eigenfunctions gye of the weakly coupled point 
scatterer. However, the result holds for any increasing sequence of numbers which 
interlaces with the norms {rij}. In particular it also applies to the eigenfunctions 
<7 A str. of a strongly coupled point scatterer. We state the theorem in full generality. 

Theorem 4.1. (Rudnick-U., 2012) 

Let T 2 = R 2 /2ttC be a general flat torus. Let a € C°°(T 2 ). Recall that g x = 
Ga/||Ga||2 denotes the L 2 -normalised Green's function. For any increasing sequence 
of numbers A = {Xj} which interlaces with the norms Af = {rij} there exists a 
density one subsequence Aoo C A such that 

(4.7) / a\g Xj \ 2 dfi -> - \^ I a dfi 



T 2 3 area(T 2 ) 



T 2 



as j — ^ oo along A c 



Remark. In particular this theorem implies that both in the weak and strong 
coupling regimes the eigenfunctions equidistribute in position space. Also note that 
the result holds for both rational and irrational lattices. The analogous result was 
proved in [26] for the standard 3D torus and general 3D tori which satisfy certain 
irrationality conditions. 

The proof of Thcorcm l4.1l uses the explicit formula (I2.15P for the Green's function 
and by approximating the Green's function by a sum over a polynomial size interval 
the problem can be translated into a number theoretical problem about the well- 
spacedness of lattice points in thin annuli. 

In a recent paper [TS] Marklof and Rudnick have shown that for rational polygons 
a full density of eigenfunctions of the Laplacian equidistributes in position space. 
In particular their result applies to the square torus. Their proof uses Egorov's 
theorem. 

Quantum ergodicity does of course not hold for the Laplacian eigenfunctions 
on the torus. If one chooses an eigenbasis of plane waves, the eigenfunctions are 
obviously localised in momentum. It is an interesting question to ask if the presence 
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of a (weakly or strongly coupled) point scatterer can change this. Curiously, the 
answer depends on the arithmetic properties of the lattice Co- 

For a weakly coupled point scatterer it is possible |12j to prove quantum ergod- 
icity for Cq = Z 2 and it is likely that one can generalise this to any rational lattice. 
On the other hand one can disprove quantum ergodicity in the irrational case. The 
failure of quantum ergodicity in the closely related case of the Seba billiard (which 
has an irrational aspect ratio) was already conjectured by Berkolaiko, Keating and 
Winn [2J. 

Open problems: It would be interesting to try to prove the analogous statement 
for the strong coupling limit: Can one disprove quantum ergodicity for irrational 
lattices? Does quantum ergodicity still hold for rational lattices? 

The microcolal lift dfij of the measures \g\f\ 2 d^L is defined by the identity 

(4.8) (Op(a)g x *,g x A= [ ad N , aeC°°(S*T 2 ). 

\ 3 3/ J S , T 2 

The quantum limits are the limit points of the sequence d\ij in the weak-* topology. 
In |10j Jakobson classified the quantum limits for the Laplacian on the square torus. 
One can pose the same problem for a point scatterer. 

Open problems: What are the quantum limits for a point scatterer on the 
square torus in the weak and strong coupling regimes? How about general rational 
and irrational 2D tori? How about 3D tori? 

5. Spectral statistics 

One of the main observations in Seba's paper |19j was level repulsion for a point 
scatterer in a rectangular billiard with Dirichlet boundary conditions and irrational 
aspect ratio. Seba computed the eigenvalues numerically and plotted the spacing 
distribution which did indeed reveal some form of level repulsion and he initially 
conjectured that the spacing distribution should coincide with that of the Gaussian 
Orthogonal Ensemble (GOE) in Random Matrix Theory. 

As Shigehara [21] discovered later, level repulsion, as observed by Seba, could 
only be expected in the strong coupling regime, where the extension parameter is 
normalised in a suitable way as the eigenvalue tends to infinity. It is likely that 
Seba carried out a truncation of the spectral equation (similar to the one in Lemma 
13. ip when performing the numerics, so effectively he calculated the eigenvalues for 
the strong coupling quantisation. 

Based on careful numerics, Shigehara predicted that one should recover Pois- 
sonian level statistics in the weak coupling regime, as was expected for the un- 
perturbed Laplacian on an irrational rectangle. In other words, the effect of the 
scatterer in this regime was too weak to have an impact in the high energy limit. In 
3D, however, the situation is very different. Numerics suggest that the eigenvalues 
of a fixed self-adjoint extension obey intermediate level statistics, just as in the 
strong coupling regime in 2D pQ. 

In a 1999 paper [4] Bogomolny, Gerland and Schmit argued that the spacing 
distribution was close to a semi-Poissonian distribution. The Seba billiard thus 
belongs to an intermediate class of systems which shows a weaker form of level 
repulsion than chaotic systems. Another example of such intermediate systems 
are flat surfaces with conic singularities. Intermediate statistics have also been 
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observed for certain families of quantum maps [8 and near the transition point in 
the Anderson model in 3D (cf. [7], section 6. 1. 3., p. 332-3). 

There are many open questions regarding the spectral statistics of point scatter- 
ers on flat tori. Even in the case of the unperturbed Laplacian, little can be said 
rigorously. Let 

Sj = rij — tij-i > 

and define 

N Co (x) = #{n j <x}. 
If Cq = Z 2 , we have the following asymptotic which is due to Landau [13] 

Br, 

(5.1) N x *{x)~^=, 5 = 0.764... 

If Co is an irrational lattice, then the multiplicity of the Laplacian eigenvalues is 
on average 4 (for a generic vector £ e £ [i] the other three choices are — £, £ and 
— £). Indeed we have in this case 

(5.2) N Co {x)~~x. 

We define the mean spacing up to a threshold x by 



71 ■ , < X 



and observe that 
(5.4) (S j 



— , if £q is irrational 



Nc (x) 



B^y/hgx, ifC = Z 2 

We introduce the mean normalised spacings 

S 

(5.5) 5j = j^, rij < x. 

Analogously let SJ = A^ 3 — l > and we introduce the mean normalised spacings 

(5.6) % = W Uj ~ X 

and it is easy to see that (Sj) ~ (^j) x as x ~ * 00 ( recan that SJ are the spacings 
of the new eigenvalues which interlace with the norms rij). 

We have the following conjecture which is due to Berry and Tabor [I] and is 
expected to hold for a generic classically integrable quantum system. We will only 
state it in the special case of flat 2D tori. 

Conjecture 5.1. Let Co be a Diophantine irrational lattice. Then the mean nor- 
malised spacings {Sj\n*<x have a Poissonian distribution of mean 1 in the limit as 
x^oo. So for any h (fC ,00 (R + ) 

(5.7) lim — V h(L) = / h(s)e- s ds 

x ^oc N Co (x) n % Jo 

Remark. We even expect this conjecture to hold in the case of the square torus 
(Important: we are ignoring multiplicities). 



QUANTUM CHAOS FOR POINT SCATTERERS ON FLAT TORI 



11 



This conjecture remains out of reach, even a proof of the existence of the spacing 
distribution is not available. However, some progress has been made with regard to 
the pair correlation function [5J 114) , in cases where the aspect ratio satisfies certain 
diophantine properties. 

It is an interesting problem to study the effect of the scatterer on the distribution 
of the spacings SJ. We introduce the difference between Laplacian eigenvalues and 
the neighbouring eigenvalue of the point scatterer: 



and we denote its mean by 



dj = rij — Aj > 



1 



(5-8) <*>. - L, 

rij <x 

One has the following bound on the mean of dj which shows that on average the 
eigenvalues rij and A J "clump" together. The result is derived in [T5] and the main 
tool is an exact trace formula for a point scatterer, similar to the trace formula 
proved in [25] . 

Theorem 5.2. (Rudnick-U., 2012) 

Let Co be a general lattice. Fix (p € (— 7r,7r). We have the bound 

(5-9) (d j ) x « 

x log a: 

Remark. This implies that if either the limiting distribution of 5j or 5j exist, then 
they must coincide. By use of the trace formula proved in [25] it is possible to extend 
this result to surfaces of constant negative curvature. 

The situation is very different in 3D. Denote £g = Z(o, 0, 0) © Z(0, b, 0) © 
Z(0, 0, 1/ab), a, b > and consider the torus T 3 = R 3 /27t£q. An orthonormal 
basis of Laplacian eigenfunctions is given by the exponentials 



M*) = 7 ^e^, ^£ 3 



where £ 3 is the dual of £§, and the eigenvalue is given by the norm |£| 2 . We denote 
the set of distinct norms by A/3. 

Let x € T 3 . In 3D the operator -A = -A|rj , where D = C£°(T 3 \ {a;o}), 
has deficiency indices (1, 1) □ and we denote the self-adjoint extensions by — A„. 
The new eigenvalues of —A v (which interlace with the norms A3) are denoted by 
77J and we define the differences 5j = nj — Tij-i and dj = nj — rn . As opposed 
to Theorem 15.21 the eigenvalues of the point scatterer on a 3D torus seem to lie 
on average in the middle of the neigbouring Laplacian eigenvalues. We have the 
following result which is consistent with intermediate statistics. The proof is given 
in [18] and the methods are the same as in the case of Theorem 15.21 



(5-10) (dj) x ~ - (*A 



Theorem 5.3. (Rudnick-U., 2012) 

Fix (p € (— 7r, 7r). For a general lattice Cq we have the asymptotic 

1 
2 



4n dimension d > 4 the analogous operator — Aq is essentially self-adjoint. 
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It would be desirable to obtain information about the variance of dj . Theorems 
15.21 and 15.31 are proved using a trace formula for the point scatterer of the type 
proved in |25j . However, this approach fails for higher moments. 

Open problems: Consider a 3D torus. Can one obtain an asymptotic for (d^J 

or a useful bound such as (dP^f < Co (Sj) x for some cq < i? What about higher 
moments? 
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